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ABSTRACT 

Results from ground-based experiments using hand-held infrared 
radiometers had shown that remotely sensed surface temperatures could be 
used to estimate surface temperatures of bare soil and to estimate crop 
water stress* An intensive ground and aircraft program was planned, in 
conjunction with the HCMM, to extend the ground-based methods to aircraft 
and satellite altitudes* The HCMM was scheduled for a one year mission 
beginning with launch in October 1977* A large agricultural field was 
selected in California (to be near the aircraft base and a wiiverslty) as 
the site for the experiment* The growing season 1977-78 was to be during 
the one year life of the HCMM* After the ground and aircraft program had 
begun, the launch date was delayed* This delay and the unexpected heavy 
rainfall made It Impossible to accomplish the original objectives In the 
detail desired* However, at three times late in the season, ground, low- 
and hlgh-altltude aircraft, and HCMM spacecraft data were obtained* 
Comparison of the multilevel acquired temperatures indicated that the 
soil moisture and yield prediction techniques developed on the basis of 
data acquired on the ground can be extended to satellite altitudes, pro- 
viding that problems of registration and accounting for atmospheric 
effects are adequately addressed* 

The plant growth and development data acquired over undulating 
terrain provided a data set that pointed out problems associated with 
Interpreting remotely sensed plant temperatures under different slopes and 
aspects* The south-facing slopes had the highest yields* The low areas 
were poorly drained and consequently had low yields* Weeds grew well in 
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the low areas and the vegetation temperatures were about the same as for 
the better sites* Yield models using stress-degree**day Inputs did not 
account for the weeds and hence overpiedlcted yields. Had the experiment 
been conducted In a dry year, the low areas would probably been the 
highest yielding areas with the south-facing slopes yielding the lowest 
because of u lack of water. 

For cloud-free days, groi'tiu- ased temperature measurements over bare 
soil were related to soil moisture content. Due to the abnormally high 
rainfall, water stress resulted from too much water, not from the lack of 
It as was expected. 

A theory was developed to explain how vapor-pressure deficit and net 
radiation affect the difference In temperature between the canopy and the 
air. This analysis showed why canopy temperatures can delineate crop 
water stress under conditions of low humidity but not under high humidity. 
This development could be used to predict when or where stress conditions 
may or may not be detectable using thermal sensors on spacecraft. 

Ground and aircraft derived temperatures were registered In space 
and time for seven flights. The spatial registration was within a 5-plxel 
grid. The relation between ground and aircraft data had an Intercept of 
near 0 and a slope of 1.07, Indicating that the aircraft temperatures were 
slightly warmer than ground temperatures at higher ambient temperatures. 
The same relation held 'd en ground temperatures from 16 sites were 
averaged and compared with aircraft temperatures Integrated to give one 
val'je for the entire 260 ha field. This result Indicates that methods 
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derived from ground-based experiments for estimating soil moisture and 
crop water stress are adaptable to aircraft programs* 

To gain an insight Into appropriate pixel sizes needed for undulating 
terrain, the 2 m x 2 ra pixels from one aircraft data set were averaged to 
produce 200 m x 200 m, 400 m x 400 m, and 800 m x 800 m pixels* When one 
high altitude data set was Included, It was found that -<lxel sizes ranging 
from 70 m x 70 m to 800 m x 800 m yielded essentially the same temperature 
Information for the entire barley field. For surveying large areas of 
grain, high resolution Is not necessary, except for the problem of the 
location of specific sites* It was extremely difficult to locate the 
experimental field, which was about 9 HCMM pixels In size* 

Multilevel temperatures, acquired from the ground, low and high alti- 
tude aircraft, and the HCMM spacecraft were compared for two day and one 
night overpasses* The U-2 and low altitude temperatures differed by less 
than 0.5*C. The HCMM data were analyzed using both the pre- and post- 
launch calibrations. The pre-launch calibration temperatures were con- 
siderably closer to the aircraft data than the post-launch calibrated 
temperatures. For the pre-launch calibrated data, HCMM temperatures were 
higher at night and lower during the day when compared to the aircraft 
data. These differences may result from the HCMM radiometer calibration 
or from the Inability to account exactly for atmospheric effects* 
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. Aerial photo of field taken on 10 March 1978. Variations In 
plant density are evident, as well as the headquarters area 
and the roads and paths. 

. The p.m. minus a.m. surface temperature difference normalized 
with the air temperature difference (top) and the surface 
minus air temperature- difference (bottom) as a function of 
gravimetric water content. The solid symbols represent dear- 
day data. 

. Theoretical relationship between the canopy-air temperature 
difference and the vapor pressure deficit. Numbers at the 
end of each line Indicate the value of the canopy resistance 
(r^) used for the calculations. Point B represents a data 
point for which a value of the crop water stress index (CWSl) 
can be obtained by ratiolng the distance BC to AC. See text 
for math.-^matical derivation of the CUSI. All calculations 
were for air temperature (T^) of 30*C, net radiation (Rn) of 
600 Wra”^, and an aerodynamic resistance (r^) of 10 sm"^. 

I. The crop water stress index (CWSI) as a function of days 

aft* r planting for site Tl. The lines were drawn by eye with 
the slope change occurring when senescence began. 
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6. Crop yie’ ' as normalized for the total receipt of solar 
rat.i?.tlon 1;-, -Ing the vegetative period )f growth versus the 
summation of stress-degree-days accumulated during the 
reproductive period of growth. 

7. The temperature difference between measurements made holding 
che Infra^'id thermometer at an angle of about 30* from horl- 
zonr ;i and when holding the Instrument vertically, as a func- 
tion of days after planting. 

3. JA.i’r versus ground (nadir) measured temperatures. 

Aircraft .lata were registered to the ground sites in time and 
space. 

9. Aircraft versus ground (nadir) measured temperatures. 

Aircraft data were averaged over entire field. Ground data 
are averages of the 16 sites. 

10. Maximum, minimum, mean, and standard deviation for 2 m x 2 m 
pixels averaged to 200 m x 200 m (4 ha) pixels (top) and 
averaged to 400 m x 400 m (16 ha) pixels f bottom). Data for 
aircraft flight at 1305, 16 May 1978. 

11. Maximum, minimum, mean, and standard deviation of 2 m x 2 ra 
pixels averaged to 800 m x 800 m (64 ha) pixels (top) and 
averaged to 1600 m x 1600 m (256 ha) pixels (bottom). Data 
for aircraft flight at 1305, 16 May 1978. 

12. Mean temperature of pixels having the highest and lowest tem- 
peratures as a function of pixel size. Lines were drawn by 
Interpolating plots of the mean temperature versus the 
logarithm of pixel side length. 
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Temperature image of the Dunigan agro-meteorologlcal experi" 
meat test field for 16 May 1978, in conjunction with a U'^2 
and HC^L1 overpass, approximately 2 m x 2 m pixel size. 
Temperature image of the field obtained at 0244, 20 May 1978, 
in conjunction with an HCMM overpass, approximately 2 m x 2 m 
pixel size. 

Temperature image of the field at 1333, 20 May 1978, in con- 
junction with an HCMM overpass, approximately 2 m x 2 m pixel 
size. 

Temperature image of the field at 1337, 16 May 1978, obtained 
with the U-2, approximately 70 m x 70 m pixel size. 
Temperature image obtained with the HCMM satellite at 
approximately 1330, '6 May 1978, approximate pixel size 600 m 
X 600 m. 

Temperature Image obtained with the HCMM satellite at 
approximately 0230, 20 May 1978, approximate pixel size 600 m 
X 600 m. 

Temperature Image obtained with the HCMM satellite at 
approximately 1330, 20 May 1978, approximate pixel size 600 m 
X 600 m. 

Satellite versus aircraft measured surface temperatures. The 
square symbol indicates the U-2 derived temperatures. The x's 
and circles represent HCMM data using the pre- and post- 
launch calibrations, respectively. 
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INTRODUCTION 

In th« early 1970*a soil physics rssssrch st ths U. S. Wstsr 
Conssrvstlon Lsboracory focused on ths diurnsl soil ttolsturs changes 
occurring In the surface layer (>10 ca) of field soils. These changes, 
driven by aoisture and temperature gradients within the soil and by 
■eterological fac:ors above the soil, required intensive field experi> 
aents to document. Although the original objectives of the experiments 
were to provide data for testing theories of simultaneous flow of heat 
and moisture in soil, they proved to be a comprehensive data base for 
testing soil moisture models used in remote sensing, and for actually 
relating temperature data that could be obtained remotely, to soil 
moisture. Some of these data were used by NASA to help show the possibi- 
lities of the HCMM satellite to estimate soil moisture. 

After developing the data to the point that it appeared feasible to 
estimate soil moisture (in the first few cm of the surface) from remotely 
sensed surface temperatures (both ground and aircraft based) we began to 
look at the possibilities of using surface temperatures to evaluate crop 
water stress and thus Indirectly measure soil moisture. While these 
ideas were in their Infancy, the opportunity to propose an HCMM experi- 
ment arose. We felt that full utilization of our work on soil moisture 
would not evolve until it was demonstrated that the methods, or modifica- 
tions of them, could be used with spacecraft data. Furthermore, with our 
emerging hopes that similar methods ccuH be developed for detecting crop 
water stress, we felt that a satellite with a thermal sensor aboard. 
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could bocoM a uaoful tool for ••••■■Ing regional crop water atatua aa 
wall aa aoll anlatura of relatively bare flelda* Thua, our HCMM propoaal 
dealt largely with eatlmatlng aoll nolature with crop t«ater atreaa being 
a minor part. Our reaearch on crop water atreaa continued at a rapid 
rate and, after our propoaal waa accepted and detailed planning for the 
experlaenta proceeded, It waa evident that crop water atreaa would be the 
major thruat of our experlmenta. Specific objectivea given In the atate~ 
ment of work (written aometlme after the propoaal was accepted) were: 

1. Evaluate the feaalblllty of monitoring aoll moisture at a selected 
bare-aoll test site on a repetitive basis by thermal methods. 

2. Evaluate the feaalblllty of monitoring plant water stress at a 
plant-canopied test site by measuring plant and air temperatures. 

3. Evaluate the feasibility of monitoring wheat biostages at a selected 
teat site by albedo and plant temperature methods. 

4. Evaluate the feasibility of estimating wheat yields at a selected 
tent site by albedo and plant temperature methods. 

'fhe long-range objective was to develop techniques for monitoring soil 
moisture, plant water stress, and crop yields on a global basis. 

The proposal was accepted in March 1976 and detailed planning for 
the ground and aircraft experiments began immediately. Launch of the 
HCMM was scheduled for October, 1977, with a planned operating life of 
one year. We were faced with the problem of conducting ground experi- 
ments over a large enough site that would be covered b>' a minimum of 9 
pixels in an HCMM image. An intensive low altitude (« 1 km) aircraft 
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program was designed to interface the ground and spacecraft data, and to 
provide data to calibrate the spacecraft sensor if necessary. The 
aircraft ptogram was made possible by VkSk/Anes Research Center. An 
excellent working relationship had developed between the NASA group and 
the U. S. Water Conservation' Laboratory during earlier cooperative 
experiments. 

We chose a site west of Sacramento, CA (38° 49'N, 121° 59'W) that 
Included about 260 ha of nonlrrigated wheat. The topography of the field 
was varied so it was expected that areas would be stressed and other 
areas would have adequate water. California was undergoing a severe 
drought during this time (1976-1977). Our plans were to begin the 
experiment in the fall of 1977 and carry it through the growing season. 
Planting was scheduled for the normal November-December time, in close 
conjunction with the scheduled October launch of the HCMM. 

In August 1977, we transferred a scientist and a technician to 
Davis, CA, for a one-year period. They made the final site selection and 
gathered the equipment necessary for the intensive ground measurements 
that were required. (Pre-dawn and afternoon HCMM time IR temperature 
measurements were to be made at 16 locations every day during the growing 
season.) About this time word reached us that the launch date had 
slipped a few weeks. Due to the severity of the drought, the farm owner 
decided to plant barley instead of wheat because of its greater drought 
tolerance. Planting began on sche'-lule. When about half of the 260 ha 
were planted, the drought broke. Rains delayed planting and caused some 
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flooding in the lower areas. Above average rainfall occurred during the 
entire season, and further delays In the launch of the HCMM were 
announced . 

We proceeded with the ground measureaents and the aircraft program. 
(The aircraft was grounded for repairs during April.) The HCMM, launched 
on 26 April 1978, took Its first data while the site was being harvested. 
Multistage data (ground, aircraft, spacecraft) were obtained for three 
satellite overpasses. 

The 670 nm of rain during the period November 1977 through April 
1978 were exceeded In this area only once since rainfall records have 
been kept. The soil was so wet that stress due to too much water was 
observed. Instead of due to lack of water as we expected. Along with the 
rain, cloudy weather prevailed much of the time, making the Interpreta- 
tion of temperature data difficult. The abnormal weather and the delayed 
launch precluded the fulfillment of several objectives. In spite of the 
problems encountered, some results of Interest were obtained. 

EXPERIMENTAL PREPARATIONS AND PROCEDURES 

Site Selection 

The criteria developed for site selection were: 

1. Terrain typical of major grain-growing areas. 

2. Good farmer cooperation. 

3. At least 250 ha area of land, nearly square, all planted to the 
same crop. 

4. Climate favorable for dryland grain production. 
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5* An area where there was a possibility of some crop wate* tress 
developing because of land relief and/or previous climatic drought 
history. 

6. Close to NASA/Ames Research Center in order to minimize flight 
times for the aircraft program. 

7. Close to a University in order to obtain essential personnel, 
space, services, and equipment. 

The location selected was a 260 ha field near Dunnigan, 35 miles 
northwest of the University of California, Davis, Yolo County, 

California. The field encompassed parts of sections 9, 15, and 16 of 
Range 1 West, Township 11 North; coordinates 38* 49*N, 121* 59’ West. 

The terrain ranged from flat to A5Z slopes. A narrow unpaved dirt road 
ran through the center of the field. 

The cropping history of the field was wheat or barley, with alter- 
nate years fallow. California had undergone a second year of drought 
and dry conditions were predicted through 1977. The owner-operators had 
previously cooperated with the University of California on other 
experiments. They agreed to per^o-mi all normal cultural operations, and 
we, in turn, agreed to reimburse them for damages caused by our experimen- 
tal procedures. 

Experimental Sites 

Sixteen sites were chosen within the 260 ha field at which the 
ground measurements would be made (figure 1). The letter R designates 
remote sites which were not fully instrumental (with the exception of R1 
near headquarters). The main fully instrumented sites were located near 
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Figure 1. Diagram of the approximately 260ha agricultural field showing 
the 16 measurement sites, roads, and paths for reaching the 
sites. 


( 
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the center of the field. Tie location symbols were; T " top of knoll, 

B ■ bottom of gully, E east-facing slope, W * west-facing slope, 

S ■ south-facing slope, and N * north facing slope. An aerial photo taken 
on 10 March 1978 Is shown in figure 2, and with the aid of figure 1, the 
16 sites can be located on the photo. The paths used to go from one site 
to another are readily detectable in the photo. Also the rolling nature 
of the terrain, the drainage patterns, and the differing plant density can 
be discerned. A description of the slope, aspect, and planting row direc- 
tion for each site is given in Table 1. The 16 sites were divided into 
two slope categories: 0-10% and >10%. Similarly, aspects (directions in 

which the slopes faced) in the four cardinal directions (north, east, 
south, and west) were selected within each slope class. 

In order to instrument several sites for intensive measurements, 
seven sites were located such that two transects (east-west and 
north-south) that intersected at right angles were formed. This was 
called the main site. 

An effort was made to space the 16 sites throughout the entire 
field, so the comparisons could be made between ground, aircraft, and 
satellite data. However, time constraints dictated that the sites be 
located so that measurements on all could be made within a 30-minute 
period. The final site selection, as shown in Figure I, came close to 
meeting these requirements. The path connecting the sites was 3.5 miles 
long and traversed slopes from 0 to 36%. Each site consisted of a 15 m 
square target area with an access path to the center from which tem- 


peratures and soil-water contents were measured 
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Figure 2. Aerial photo of the field taken on 10 March 1978. Variations In plant 
density are evident, as well as the headquarters area, roads, and paths. 
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Table 1* Identification symbol » slope, aspect, and direction of planting 
for the 16 ground measurement sites* 


Slope 

(X) 


Aspect 

(Deg.) 


General Planting 
Direction 


Prlman 


Overplantlni 



R1 

2 

342 

NE-SW 



R2 

22 

52 

NW-SE 



R3 

24 

247 

N-S 



R4 

20 

157 

NW-SE 


i 

T1 

0 

- 

E-W 

NW-SE 


W1 

4 

232 

E-W 

NW-SE 


B1 

3 

322 

N-S 


* 

El 

7 

62 

N-S 



SI 

3 

162 

N-S 

NW-SE 


B2 

2 

160 

ME-SW 



N1 

b 

347 

NE-SW 



R5 

36 

47 

NW-SE 



R6 

24 

7 

N-S 


r 

R7 

0 

- 

E-W 



R8 

0 

- 

E-W 



R9 

1 

132 

N-S 
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The soil* are of the Sehorn-Balcom complex* They are well drained, 
gently to ateeply eloping silty clay loams and clays on sandstone* The 
Balcom soil Is classified as a fine-loamy, mixed, calcareous, thermic, 
typlc Xerorthent, while the Sehorn Is a fine, montmorlllonltlc, thermic, 
entlc Chromoxerert * Soil samples taken prior to planting showed a 
variable but relatively low amount of nitrogen In the upper 25 cm of 
soil* 

Facilities 

Since the field was some distance from all services (power, etc*), 
the experiment had to be completely self-contained* Primary power was 
furnished by a 32 KW diesel generator* A portable building and two 
mobile trailers (which housed the electronic data logging and teat 
equipment, oven, refrigerator, balance, battery chargers, and other 
miscellaneous devices) received power from the generator* To Insure con- 
tinuity in data collection, th*' data logger obtaii.^d power from storage 
batteries when the generator failed* The data logger could operate for 
2-3 days without the generator* A 15 KW propane generator was available 
as a backup for the diesel generator* Diesel and propane were delivered 
by local fuel companies as needed* A radio telephone was Installed for 
emergency purposes* 

Measurements 

Ground-based measurements %rare divided into two types: automatic 

and manual* Automatic refera to measuring several micrometeorologica*. 
parameters near the trailers (Rl) and at main site (Tl, Wl, Bl, El, 

SI, B2, and Nl) and recording the 20-mlnute averages of these parameters 


r 
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for tho tntlro duration of tha axparlaant. Theaa paraaatara and tha altaa 
at which chay wara aaaaurad ara glvan in Tabla 2. 

Tha manual naaauramanta ara thoaa which could not ba racordad on a 
data logger, and %fare taken on a dally or leas frequent basis* Table 3 
describes the parameters measured, the locations, and the frequency of 
measurement* Table 4 gives information on the instruments used* 

Crop and/or soil surface temperatures were measured with a portable 
Infrared thermometer by going to the center of the plot and aiming the 
instrument at a spot about 3-5 m from the center in each of the four car- 
dinal directions* Also, a down-looking (nadir) temperature was taken at 
four points at the plot center* In this way eight temperatures were 
obtained at each site, four at an angle and four at nadir* 

Temperatures were measured every day from 17 October 1977 through 
23 May 1978 except in inclement weather or if an equipment malfunction 

I 

occurred. The data were taken about one hour prior to sunrise and again ^ 

1 

about one hour after solar noon* This amounted to 256 Individual 

infrared thermometer readings daily* Addltonally, dry bulb and tiet bulb I 

! 

air temperature measurements were taken at each of the 16 sites during \ 

1 

each of Che twire-a-day measurement periods* I 

Soil tiater contents were meartred with t neutron scattering probe 

and scaler, three times each week* Measurements %«re made In 20-cm 

* \ 

increments to a maximum depth of 160 cm. In cases where the soil was 

J 

shallow or a water table was present, the maximum measurement depth was 

less chan 160 cm* I ; 

i \ 
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Table 2* Location of Ipatrumanta uaad for automatic data logging of 
metarologlcal paraaatara* 


Parameter 


R1 T1 W1 


Incoming aolar radiation 


X X 


Reflected aolar radiation 


XXX 


Net radiation 


XXX 


wind speed and direction 


X X 


Air temperature @ 150 cm 


XXX 


Soil temperature 9 aurfbce 


XXX 


Soil temperature 9 2 cm 


XXX 


Soil temperature 9 4 cm 


XXX 


f 
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Frequency of oenuelly obteined ground-beeed aeeeurenente* 


Teble 3. 


Pereneter 

Surfece tenpet -\ture 
w/lnfrered rediometer, end 
dry end «ec bulb teapereturee 
with en eeplreted peychroaeter 
preeur.*lae and 9 1230-1300 hr 

Soil aoleture w/neutron 
aeter every 20-ca depth to 
160 oi 

Ten plant eeaplea froa 
each alte for height, dry 
weight, green leaf area, 
end growth at age 

Soil teaperature 
w/theraietor probe 9 2 - 
end 4-ca depth at eltee 
R2 through R9 

Albedo at sites R2 through 
R9 w/portable unit on 
clear days 

Insect sweeps and pathogen 
Identification 

Froa weather station 9 R1 
site; aaxlaua-alnlaua 
teaperatures and evapora- 
tion pan (also, weekly 
record froa hygrotheraograph)* 

Ralngages at all 16 sites 
(as appropriate) 


Sun Mon Tue Wed Thu Frl Sat 


X X X X X X X 


XXX 


X X 


X X X X X X X 


X 


X X X X X X X 




NASA aircraft at 1>4 ka to aeasure 
surface teaperature 9 presunrise 
and 9 1300 hours 


X 


X 
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Tabl« 4. Manufaetursr and 

modal of ground-baaad 

inatruments used* 


Paraaccar 

Instrument 

Manufacturer 

Modal 

Incoalng solar radiation 

Pyranoneter 

Ipplay 

8*48 

M H M 

•s 

Lambda Inst. Co* 

7407 

•• M H 


Spectran 

4048 

Raflactad solar radiation 

•s 

Bppley 

8-48 

M M M 

•s 

Spectran 

4048 

Nat radiation 

Hlniatura Net 
Radiometer 

Micromat Inst* 

— 

Wind spaad & direction 

Propeller Vane 

Gill 

35003 

Air & aril teaperature 
Manual Measurements 

Oopper-constantan 

thermocouples 

Thermo Electric 


Surface temper* rure 

Infrared thermometer 

Telateap 

44 and 
AG-42 

Dry and %#et bulb temp* 

Aspirated psychro- 
meter 

Bendix 

566 

Soil %mter content 

Neutron Moisture 
Meter 

Troxler 

2601/1255 

Soil temperature 

Thermistor probe 

Extech 

1200-1108 

Albedo 

Pyranoaeter 

Spectran 

8-48 

Maxlaua>ainimua air temp* 

Mercury Thermometers 


— 

Air temp. 4 humidity 

Hygrothermograph 

Belfort 

— 

Rainfall 

Ralngauge 

Tru-Chek 

— 

Evaporation 

Class A pan 

— 

— 


IS 


T«n plant randoaly aalactad froa aach of eha 16 altaa, aara 

tatharad twlca waakly and cranaportad to Davit for analyala* Xnforaatlon 
gatharad froa thaaa aaaplaa waa phonological growth ataga, plant haight, 
dry waight, groan loaf area, nuabor of groan laavaa par plant, nuabar of 
tilloro, and nuabar of haada* Shortly aftar aaarganca a plant danaity 
count waa aado at each aita> At thraa locationa at tha aain aita (Tl, VI, 
and SI) thara wara parta of tha target araa that wara ovarplantad (tha 
plantar want ovar it twica)* 

At tha and of tha BaaauratBant pariod, a portion of aach aita waa 
hand harvaatad. Fiva conaacutiva ‘‘aquaraa" wara takan in aach of tha four 
cardinal diroctiona, starting one natar froa ifhara we noraally stood to 
take plant tamparaturas* All plants wara ranovad froa these 20, one natar 
squares at each location, tha nuabar of heads in aach aquare aetar 
counted, and the samplas labeled and bagged for further processing. 

Following harvest, all aaaplaa ware thrashed, tha nuabar of saads 
counted, and their weight and water content detarainad. Froa these data, 
tha following inforaation was obtained: Seed water content; nuabar of 

heads par square aetar; nuabar of seeds par head; nuabar of seeds par 
square aetar; bushel waight; and yield. Tha last two itaas were 
corrected to a graviaatric seed water content of 0.125 to allow coa- 
parisons between sites. 

Albedo, tha ratio of reflected to incoaing solar radiation, was 
aaasurad with a pyranoaatar attached to a pole %rith tha sensor facing tha 
crop. Tha davJca was hand-held ovar tha crop and asasureaants were aade 
at six locations around the footpath leading to tha canter of tha raaota 
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f ploec* locoBlnf solar rsdiatioa was Mssursd cwlcs at aach plot* Tha 

■aasuroMnea wars aada naar solar noon on elaar days* 

In ordar to asaass crop damaga dua to inaacts, swaap nat saaplaa 
( vara cakan ones a waak and insacts wars countad and Idanciflad* Tha top 

20-30 CB of Cha plants wars saaplad using 25 swaapa with tha nat* 
Approxinataly 12*5 wara saaplad at aach sita, onca aach waak. 

^ Pathogan activity was notad and idantifieation was made by tha county 

agricultural agant* 

A routine waathar station was oaintainad northeast of tha haad- 
quartars trailar* Daily obsarvationa of oaxiaua and ainiaua air taa- 
paraturaa and pan avaporation wara racordad* Alao» waakly racorda from a 
hygrotharoograph wara collaetad* 

Raingagas at aach of tha 16 sitaa wara read after aach rain* 

Field cultural operations 

Tha field had bean diac-plovad during tha sunoar 1977. Tha faroars, 
anticipating another year of drought, planted Briggs barley (Hordaua 
vulgara i**) instead of wheat (tha drought tolerance of barley is greater 
than wheat), during tha week of 1-7 Dacaabar 1977 at tha rata Oi 112 

/ 

kg/ha. They also incorporated 35.6 kg/ha of nitrogen and 39.4 kg/ha of 
phosphorous fertiliser in tha fora of aaaoniua phosphate* On 1 February 
1978, 44.8 kg/ha additional nitrogen was applied by aircraft, because ouch 

t 

of tha original for' .liter had bean leached froa tha root tone by heavy 
rains. During the following 2 weeks, tha barley turned dark green. Areas 
alssad in the 1 February application wara aerially fertilised with 112 
kg/ha nitrogen in tha fora of aaaoniua nitrate on 18 February. On 
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20 Febmary 1978» a herbicide was aerially sprayed to control broadleaf 
weeds* From 23 to 30 May 1978, the barley was harvested. 

Weather 

The weather during the growing season was characterized by above 
normal rainfall and temperature. A summary of the rainfall Is given in 
Table 5 for each site by months. During October and November 1977 there 
were only two ralngages Installed, so the numbers in parentheses are 
estimated based on the relationships between gages for later months. It 
turned out that this year Indeed marked the end of the drought that 
California had experienced for several years. In fact, our experimental 
site received nearly 70Z more rain (675 nmi) than normal (400 mm) during 
the 7-month period of measurement. 

The overabundance of rain, the second wettest year on record, caused 
several problems. First, crop water-stress conditions did not develop at 
any of the sites. Second, the large number of Inclement days resulted in 
fewer aircraft overflights than were planned. Third, the excess water 
caused severe flooding and erosion problems* Five of the sites, Rl, Bl, 
B2, R7, and R9 developed water tables that were observable In the neutron 
access tubes. Erosion in the gullies caused two sites, B2 and R7, to be 
relocated, creating some discontinuity in the data. All sites with high 
water tables had poor stands of barley, and a severe weed problem. 

Low Altitude Imagery 

Low altitude thermal imagery of the site was acquired from a Cessna 
402 aircraft. Data were taken throughout the growing season, from 
planting to harvest, except during April when the aircraft was grounded 
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Table 5> Summary of rainfall at the sites for the period October 1977 
through ^rll 1978. 


Oct 77 


Gate 

4.6 

HQ 

(5.6) 

R1 

5.6 

R2 

(5.6) 

R3 

(5.8) 

R4 

(4.6) 

T1 

(5.6) 

U1 

(5.3) 

B1 

(5.3) 

El 

(5.3) 

SI 

(5.3) 

B2 

(5.6) 

N1 

(5.1) 

R5 

(6.1) 

R6 

(6.6) 

R7 

(5.3) 

R8 

(5.6) 

R9 

(5.3) 

X 

5.6 

Zx 

5.6 


Nov 77 

Dec 77 

(mm) 
Jan 78 

76.5 

70.6 

229.4 

(86.9) 

77.5 

262.9 

80.3 

67.8 

230.6 

(85.6) 

74.2 

254.0 

(88.4) 

74.2 

267.2 

(68.6) 

57.1 

210.1 

(87.4) 

75.9 

262.1 

(82.0) 

70.1 

247.1 

(82.3) 

72.4 

245.9 

(83.8) 

69.6 

248.9 

(82.0) 

73.4 

242.1 

(85.3) 

71.6 

255.0 

(78.7) 

67.6 

234.2 

(92.2) 

82.6 

272.3 

(100.3) 

82.8 

301.5 

(82.6) 

69.9 

246.6 

(87.6) 

74.7 

264.2 

(81.8) 

71.9 

241.3 

84.1 

72.4 

251.0 

89.7 

162.1 

413.0 


Feb 78 

Mar 78 

Apr 78 

111.5 

84.3 

54.1 

116.6 

84.9 

54.1 

113.5 

84.1 

55.6 

122.2 

88.6 

56.9 

123.4 

93.2 

58.7 

93.7 

77.5 

53.8 

120.9 

83.1 

58.9 

113.5 

82.8 

57.1 

114.3 

84.1 

57.4 

121.9 

83.8 

56.4 

115.6 

83.6 

56.6 

121.9 

85.6 

56.4 

112.5 

83.6 

58.9 

129.3 

93.5 

59.7 

142.2 

94.5 

63.0 


93.2 


116.8 

87.1 

58.7 

120.9 

86.1 

56.6 

116.3 

86.3 

54.4 

118.1 

86.4 

57.2 

531.1 

617.5 

674.6 
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for maintalnance. Hore than 60 flights were completed, of which 55 pro- 
vided usable data. Flights were made both prior to sunup and about one 
hour after solar noon. Thermal Imagery was acquired with a Texas 
Instruments Model RS-25 Infrared scanner operating in the 10.5-12.5 wm 
bandpass region. The Instrument had an Instantaneous field of view 
'IFOV) of about 2 x 2 m at a flight altitude of about 1.2 km (4000 ft.) 
and a temperature accuracy of about 0.2"C. The scanner contained two 
blackbody calibration sources with platinum resistance thermometers for 
continuous Inflight calibration. At the completion of each flight, 
atmospheric temperature and humidity were measured at various levels down 
to near ground level. These were used to correct the thermal-IR data for 
absorption by atmospheric water vapor. 

The RS-25 scanner data were recorded on magnetic tape in analog 
form, converted to digital form, and processed by using digital image 
techniques. The first step in converting the analog data to a computer- 
compatible digital format was to average it using "sample and hold" 
techniques, which Integrate the analog signal for a predetermined period 
of time. The Integrated signal was then measured by an A/D converter 
which assigned a value to represent the signal's relative amplitude. The 
"sampling" period represents one picture element (or pixel), and the 
pixel amplitude is converted to a digital value ranging from 0 to 255 (8 
bits). The thermal Infrared data was processed on an HP-3000 computer 
system equipped with a COMTAL video display, two 1600 BPI tape drives, 
and two 50 megabyte disk memories. The software program, called IDIMS, 
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acted upon the digital Image* Since the RS-2S scanner collected data In 
a scan-llne-to-llne sequence, each line represented a row In the digital 
Image. Similarly, the selected sample Interval for digitizing determined 
the number of pixels In each scan line. Since each scan line contained 
the same number of picture elements, the pixels made up the columns In 
the Image array. 

Seven separate processing steps were used to generate the final 
output images reported here. They were: 

(1) Scene selection. All recorded digital Image data were 
reviewed on the COMTAL display. Thermal calibration, high frequency 
Image Jitter, bad or missing scan lines, image noise, and general Image 
appearance were used as selection factors* 

(2) Scene reduction. The portion of the scene outside the area of 
Interest was edited out. This conserved processing time. 

(3) Sweep distortion correction* Geometrical distortions caused 
by the constant angular velocity of the scanner mirror were corrected* 

(4) Transfer of black body digital values to the reduced scene. 

This step was required since the position of the calibration signals on 
the original total scene caused them to be edited out In Step 2, above. 

(5) Thermal calibration. The basic digitized Image was created 
with 256 (8 bit) grey levels. These levels were proportional to the 
energy received by the RS-25 scanner from each point on the ground* The 
thermal black body references were also imaged within 256 grey levels, 
but their temperatures were known* Thus, all Image grey levels could be 
transformed Into apparent temperatures by using the black bodies as 
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function generators* The output image, then, had real-valued pixels 
which represented the apparent temperatures on the ground as seen by the 
RS-25 scanner. 

(6) Atmospheric correction. The foregoing apparent temperatures 
were corrected for atmospheric effects by use of the LOWTRAN 4 computer 
code. Multi-level, aircraft acquired, values of atmospheric temperature 
and humidity were utilized. The effects of CO2, O3 and aerosols were 
accounted for in the code. 

(7) Presentation of the processed Image* Final results %iere pre- 
sented In the following forms: 

a) Pseudo-colored video display, where each color represents 
a discrete temperature Interval. 

b) Total scene display or selected areas expanded to fill 
video screen size. 

c) Line printer output of apparent temperature values or 
temperature differences* 

d) Single pixel values selected by the operator. 

High Altitude Imagery 

High altitude thermal imagery of the site was acquired with the 
HCMR prototype scanner flown on NASA's U-2 aircraft at an altitude of 20 
km (pixel size 70 x 70 m). Data were acquired the afternoon of 16 Hay 
1978, simultaneous with the HCMM spacecraft pass. The data tapes from 
these flights were made computer-compatible by NASA/GSFC, and returned to 
NASA/Ames where they were processed In a manner similar to that described 
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In the previous ssctlon* Air tsnpsrsture, and dew point were acquired 
from the Oakland, California National Weather Service soundings for use In 
computing atmospheric corrections to the thermal data by the LOWTRAN 4 
computer code. 

Spacecraft Imagery 

Multi-staged experiments were conducted on the days and from the 
platform given In Table 6. 


Table 6. Dates and platforms for acquiring multi leval aircraft 
and spacecraft data over the experimental field. 


Date 

Ground 

Cessna 

W-2 

HCMM 

5/16/77 

Postnoon 

X 

X 

X 

X 

5/20/77 

Presunrise 

X 

X 


X 

5/20/77 

Postnoon 

X 

X 


X 


The spacecraft data were handled in a manner similar to the aircraft 
data except that NASA/GSFC provided computer-compatible tapes. The 
mapping of the grey levels to corresponding temperatures was done using 
a set of parameters and an equation supplied by Goddard: 

T(I) - K2/lln{Kl/(l-K3)+l]| (1) 

where T Is the Kelvin temperature, 1 the pixel byte level, Ki - 14421.587, 
K2 - 1251.1591, and K3 - -118.21378. 
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RESULTS AND DISCUSSION 

Soil Moisture Estimation 

The objective of monitoring soil moisture In the surface layer of 
bare soli on a repetitive basis was addressed In the Fall of 1977 prior 
to emergence of the barley. Soil samples were taken from the surface to 
4->cffl at about 1100 (Jackson et al.» 1976) for 16 days over a 30-day 
period. Soil moisture was determined on a weight basis and expressed as 
gravimetric water content. TWo temperature parameters were used: the 
surface minus air temperature In the afternoon (Tg - ^A^pm 
1975; Reglnato et al., 1976), and the afternoon-morning surface tem- 
perature difference as normalized by the method of Idso et al., 1976. 

The results are shown In Figure 3. In theory, both temperature parameters 
should decrease with Increasing water content. With some Imagination, one 
can see a decrease with Increasing water content for the afternoon 
surface-air temperature difference, but not for the day-night differences. 
The scatter Is the result of the adverse weather conditions. Of the 16 
days of measurement, only two were relatively clear. These two days are 
Indicated by the solid circles In Figure 3. Considering only the two 
clear days, the temperature parameters did decrease with Increasing 
water content. The two points compare well with published data (Idso et 
al., 1976). 

The scatter l ..he data In Figure 3 emphasizes the fact that ther- 
mal Infrared techniques for estimating soil water content yield reliable 
results only under clear sky conditions. Thermal Infrared data from 


ORIGINAL PAGE IS 
OF POOR QUALITY 


GRAVIMETRIC WATER CONTENT(g/g) 


Figure 3. The p.m. minus a.m. surface temperature difference normalized 

with the air temperature difference (top) and the surface minus 
air temperature difference (bottom) as a function of gravimetric 
water content. The solid symbols represent clear day data. 
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•Atallltes can only b« obtained under clear aklea* Thua, this technique 
la moat uaeful In araaa with few clouds. Such areaa normally are arid 
and, therefore, aoll molature Information la of particular Importance. 
Plant Water Streaa 

Data relating to our objective of monitoring plant water atreaa 
were collected from January through May of 1978. In the courae of our 
Inveatlgatlona It became evident that the atreaa-degree-^ay concept aa 
propoaed by Idao et al., 1977 and Jackaon at al., 1977 was an Inauf- 
flclent Indicator of stress, especially under humid conditions. 
Consequently we derived a crop water stress Index (CUSI) from energy 
balance ccusiderations (Ji^ckson et al., 1981, in press) and empirically 
from a large data set (Idso et al., 1981, in press). A summary of the 
development of this Index will be presented here not only for the pur- 
pose of quantifying plant water stress but also to provide a basis for 
Interpreting spacecraft thermal data obtained over different climate 


regions. 

The energy balance for a crop canopy can be written 

- G + H + XE (2) 

where Rf) Is the net radiation (W m~^), G is the heat flux below the 
canopy (W m~^), H Is the sensible heat flux (W m~^) from the canopy to 
the air, AE Is the latent heat flux to the air (W m**^), with X being 
the heat of vaporisation. In their simplest forms, H and E can be 
expressed as 


H " pCp (Tg “ 

XE - PCp (ej - «A)/tY(r, + r<.)J 
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wh«r« p It tha dtntlty of air (kg Cp tha haat capacity of air 

(J kg'^C^), tha turfaca tamparatura (*C), tha air taaparatura 
. a 

(*C), a^ la tha aaturatad vapor praaaura (Pa) at a^ tha vapor 
prettura of tha air (Pa), y It tha ptychroaatrlc conatant (PaC~^), r^ 
la tha aarodynamic raaiatance (a m**!-), and r^ la tha canopy raalatanca 
(a a~^) to vapor tranaport* A datallad dlacuaalon of procaduraa 
loading to aquationa (2), (3), and (4) la glvan by Hontolth (1973). 
Conblnlng (2), (3), and (4), ataumlng that G la nagllglbla, and 

daflnlng 4 at the alope of tha aaturatad vapor pieaaure ~ tamparatura 

* * . 1 
relation (e^ - a^)/(Te ~ T^), unlta of PaC*^, tie obtain 


Tc - Ta - 


taRn ’^c/’^a) 


* 

«A " «A 


PCf 


A + Y(1 + Tc/'a) A + Y(1 + fc/^a) 


(5) 


which relatea the difference between the canopy and the air temperaturea 

* 

to the vapor preaaure deficit of the air (e^ “ e^), the net radiation, 
and the aerodynamic and crop reaiatancea. 

The upper limit of can be found from (5) by allowing the 

crop realatance r^ to increaae without bound, l.e., aa r^ '»' • 


Tc - Ta 


ra*n 

PCp 


( 6 ) 


The lower bound, found by aettlng r^ ■ 0 In (5) (the caae of wet 
planta acting aa a free water aurface), la 


I 
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Tc - Ta - 


(•a “ «a) 


A + Y 


A + Y 


Equations (5) and (7) daacrlbe a linear relation between 


and the vapor preaaure deficit, e^ “ e^* Thus, for a particular 


teaperature, the lower bound ia a line extending froa the intercept at 


* * * 

*A ^ *A " ^ (eaturated air) to a value of «a ~ *A * *A (coapletely dry 


air). Since A appears in both the slope and the Intercept, both teras 


are teaperature dependent. 


Equation (7) represents the case of evaporation froa a free water 


surface, which is not necessarily the case for potential evaporation 


from a crop. In irrigated areas, the soil aay be adequately supplied 


with water, with the plant surfaces being dry. In this case the canopy 


resistance is probably not zero (van Bevel and Airier, 1968), but has a 


value that we will call the canopy resistance at potential evapotrans-* 


piration (r..). The value of r>. will probably be different for dif- 


ferent crops, and may change with crop variety. Setting r^ - r^p in (5) 


we have 


Tc - Ta 


* 

«A “ *A 


A + Y* A + Y* 


where 


Y* - y( 1 + rcp/>^a) 


A crop with adequate water will transpire at the potential rate for 


that crop. As water becomes limiting, the actual evapotranspiration 
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will fall balov tha potantlal rata* A naaaura of tha ratio of actual to 
potantial evapotranspltatlon should, tharafora, ba an Indax of crop 
vatar status* Conblning (2), (3), and (4) and solving for XB yialds 


XE - 


A Rn + PCp («A - 


( 10 ) 


4 + Y (1 + r^/r,) 

which is tha Panaan**Montaith aquation for avapotranspiratlon in tarws of 
canopy and aarodynaalc rasistancas (Montaith, 1973: Thoa and Olivar, 
1977)* Taking tha ratio of actual (XE for any r^) to potential (XEp for 

■ 'cp) 

A + Y* 


E/E, 


( 10 ) 


'P 4 + Yd + rc/r,) 
with Y^ daflnad by (9). Jensen (1974) and Howell at al (1979) discussed 
(11) for the case of r^p - 0 , l.e., Y* ■ Y* Rearranging (11) will give 
in terms of E/Ep, a result reported by van Bevel (1967), Szelcz and 
Long (1969), and Aussall (1979), again with r^p ■ 0* Van Bavel measured 
E with lysimeters and calculated the canopy resistance* 

The ratio E/Ep ranges from 1 (ample water, r^ ■ r^p) to 0 (no 
available water, r^ -^^ ■*)* In studying plant-water relations one thinks 
of a plant as going from a no-stress to a stressed condition* 

^erefore, it is esthetlcally pleasing for a stress index to go from 0 
to 1* He consequently define a crop water stress index (G^SI) as 

Y(l+ r^/r,) - Y* 


CWSI 


1 - E/Eo - 


P 4 + Yd + r^/r.) 


( 12 ) 




29 

To coleulato tho CUSl or E/Ep uolng (11) or (12) roquirct a valua 
for cha ratio r^/r^* Thla la obcalnad by raarranglng (5) vlth cha 


raault 

If fa^>/(PCp) ~ (Tc - Ta)(A » T) - («A ~ •A> 
^a ^ t(^c ” 


(Vi) 


giving tha ratio r^/r^ In tama of not radiation, canopy and air 

j 

taaparaturaa, vapor praaaura daflclt, and tha aarodynaalc raalatanea* 

In praetlca, r^/r^ la avaluatad ualng (13) and aubatltutad Into (12) 
to obtain tha QfSI* 

Tha alopa of tha saturatad vapor praaaura “ tanparatura ralatlon 
(A) appaara In aoat of tha aquatlona In tha pravloua aactlon* Aa a 
flrat approximation, A can ba avaluatad at tha air taaparatura (Tj^). 

VHian tha taaparatura diffaranca (T^ - T^) la larga (for tha caaa of 
wall-watarad cropa at high vapor praaaura deflclta) a battar approxima- 
tion la to avaluata A at (T^ Obvloualy, when T^ la near tha 

two approachaa yield almllar reaulta. Taking tha average of tha canopy 
and air tamparaturaa la aufflclant for (5), (7), (8), and (13), but (11) 
and (12) poae a problem* Following cloaaly tha development of (11) from 
(10), wa find that A In tha numerator ahould ba avaluatad at tha average 
of cha air tamparatura and cha canopy Camparature chat would obtain If 
f tha crop waa evaporating at potential . Tha A In tha denominator ahould 

I 

ba evaluated at Cha average of tha actual meaaurad canopy temperature 
and tha air tamparatura* Keeping A aa cha value at tha meaaurad tam- 
f paraturaa and A* aa Cha value ualng the calculated canopy tamparatura at 
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f r 

i 

j 



potential, the nthaorutor of (11) bacoaaa A* ^ T* and thr. noaarator of 
(12) bacoMs y( 1 + r^/r,) - y* + (A - A*). 

Tha avaluatlon of A* la eoaplieatad by tha fact that at potan- 
tial nay not ba known* Thla can ba calculated uaing an itaratlva proca** 
dura with (8) by evaluating A at calculating T^, evaluating a new 
A at (T^ •¥ Tj^)/2, and recalculating until an acceptable value la 

obtained. In practice tha uaa of (11) and (12) with A evaluated at tha 
average of tha aaeourad canopy taaperatura and the air teaparatura will 
yield alallar raaulta lor both tha low and high valuea of tha Indlcaa, 
with tha aaxlaua error occurring near 0*5. At tha aldpolnt, tha dlf* 
faranca between tha two aethoda of calculation will ba wl Mn 0*06, 
a*a*, 0.48 for tha one value of A and 0.54 whan A and A* are uaad. Tha 
raaulta reported hare are In taraa of A evaluated at both tha canopy 
teaparatura at potential and at tha average taaperatura of the aeaaured 
canopy and air teaperaturea. 

A graphic example of the CUSl la given in Figure 4. Ihe four llnaa 
eaanallng from the ordinate and ending at about 4.24 kPa were calculatad 
uaing (5) with ■ 600 Ua~^, r^ “ 10aa~^ and r^ taking on the valuea 
ahown at the end of the llnea (l.e., 5, 50, 500, • canopy rc4 atancea). 
The lower line (r^ “ 5) repreaenta a wall-watered crop tranaplrlng at 
tha potential rata. Tha topmoat line (r^ ■ *) repreaenta tha caae of a 
nontranaplrlng crop. If a canopy temperature la maaaured at a vapor 
praaaure deficit of 3 and found to be 3*C below the air temperature, tha 
datum would plot at B. If tha planta were not atraaaed the point would 
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r * Figure 4. Theoretical relationship between the canopy-air temperature 

> ' difference and the vapor pressure deficit* Numbers at the 

^ I end of each line indicate the value of the canopy resistance 

I used for the calculations* Point B represents a data 

, ^ point for which a value of the crop water stress index (CUSI) 

i ' can be obtained by ratloing the distance BC to AC* See text 

for mathematical derivation of the CUSI* All calculations 
! were for air temperature (T/^) of 30*C, net radiation (R^) of 

I 600 Um~’^, and an aerodynamic resistance (r^) of 10 sm~^* 
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be at C» and if no transpiration could take place, the point would be at 
A. The CWSI is BC/AC. 

The CWSI was calculated for all sites having valid net radiation 
data, beginning about 15 February. Data for clear days on pi''* T1 are 
shown in Figure 5. Except for two data points (days 100 and 112) the 
results show a slight Increase in the CWSI until about day 140. At this 
time the plants had begun to senesce. As senescence progressed and the 
plants matured the OfSI Increased. Data for other plots were similar to 
those shorn in Figure 5. The value of the CWSI depends on the choice of 
the canopy resistance at potential evaportranspiration (r^p). We used a 
value derived from irrigated wheat data (5 sm~^), since a value has not 
been established for barley. A larger value of Tf.p would have resulted in 
smaller values of the CWSI. Since more than adequate moisture was 
available most of Che season, and since all plots for which Che CWSI was 
calculated had similar values of the (MSI, it appears that our choice of 
r,.p was low. 

Equation (5) and Figure 4 provide an aid for interpreting thermal 
infrared data over different climatic regions. In humid regions, the 
vapor pressure deficit is low and plant temperatures are generally 
warmer than air temperatures. The range of canopy-air temperature dif- 
ferences are small and therefore thermal discrimination of different 
vegetation types or stages is difficult, if not Impossible. In an arid 
area, where the vapor pressure deficit is usually high, discrimination 
between well watered and stressed vegetation is possible. Also 


DAYS AFTER PLANTING 


Figure 5. The crop water stress index (CWSI) as a function of days after 

planting for site Tl. The lines were drawn by eye with the slope 
change occurring when senescence began. 





senesclng vegetation could be discriminated from actively growing 
vegetation. 

Plant Growth and Yield 

Plant samples were taken twice each week over the growing season. 

All plant measurements were made by personnel at the University of 
California, Davis. Details of these measurements are given In an 
excellent thesis written by Carol Whitman (Whitman, 1981) and submitted 
toward fulfillment of requirements for a Master of Science degree. 
Pertinent data abstracted here are the ttl'lr.al plant growth stages 
(Table 7), the final 3 rleld of grain and total dry matter (Table 8), and a 
composite of sites according to drainage and Irradlance (Table 9). The 
growth stage data presented In Table 7 show that the time required to 
reach heading and anthesls was different for the different sites; also 
yields were different for the various sites. For the moat part, 
differences were due to slope and aspect of the terrain, and the fact that 
several sites had high water tables during much of the season. 

Temperature measurements did not distinguish these differences. 

The crop water stress Index, because of the recentness of Its deve- 
lopment, has not been Incorporated Into yield models. The predecessor of 
Che CWSI, the stress-degree-day (SDD), has been used In ylelu models. 
Although some shortcomings of the SDD were discussed earlier In this 
report, we used the model and data of Idso et al. (1980) to show that the 
yields obtained In this experiment fall well within the range of yields 
obtained for a number of crops (figure 6). The encircled open diamonds 
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Table 7. Critical development stages, In days after planting, for barley 
grown at Dunnigan, 1977>78« Day 0 tias 7 December 1977. 


Site 

Planting 

Begin 

Tillering 

Rapid 

Tillering 

Jointing 

Boot 

Heading 

Anthesls 

Maturity 

R1 

0 

32 

54 

98 

123 

139 

146 

158 

R2 

0 

32 

51 

81 

108 

126 

134 

151 

R3 

0 

32 

51 

78 

104 

123 

132 

151 

R4 

0 

32 

51 

74 

100 

119 

127 

148 

R5 

0 

32 

51 

80 

106 

125 

134 

153 

R6 

0 

32 

51 

80 

106 

125 

134 

150 

R7 

0 

32 

51 

83 

108 

127 

135 

158 

R8 

0 

32 

51 

82 

108 

126 

134 

150 

R9 

0 

32 

51 

85 

113 

135 

144 

154 

T1 

0 

32 

51 

78 

104 

123 

131 

155 

W1 

0 

32 

51 

78 

104 

123 

131 

149 

B1 

0 

32 

51 

86 

114 

133 

142 

152 

El 

0 

32 

51 

84 

111 

129 

137 

158 

SI 

0 

32 

51 

80 

105 

124 

132 

156 

B2 

0 

32 

51 

80 

105 

124 

132 

160 

NI 

0 

32 

51 

86 

114 

130 

138 

159 
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Final yield and total dry matter for barley grown at 
Dunnlgan, Callfornia» 1977-78* 



Yield 

Total 
Dry Matter 

Site 

(gn"^) 

(gm"2) 

R1 

104.5 

364 

R2 

321.5 

1057 

k3 

372.7 

1621 

R4 

417.5 

1373 

R5 

315.3 

1049 

R6 

234.0 

883 

R7 

232.5 

599 

R8 

367.2 

1249 

R9 

118.8 

634 

T1 

325.1 

1111 

W1 

317.2 

988 

BI 

116.2 

433 

El 

279.1 

546 

SI 

348.5 

782 

B2 

— 

570 

N1 

243.5 

706 


4 

i 

I 
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Indicate data for 11 of the 16 altea. Of the five altea not shown, one 
was not hand -harvested before the fanner operated combine arrived* Four 
ocher sites were In low poorly drained areas* These sites had a poor 
stand of barley but a lush stand of weeds* The foliage (In contrast to 
crop) temperatures were low, consequently the SDD's were low* The SDD 
yield model would predict high yields under these weedy conditions, 
whereas the actual barley yields were low* These results demonstrate that 
weed growth can be a serious confusing factor when canopy-temperature- 
based yield models are used* 

Although the yields from 11 sites 'fall within the range shown In 
figure 6, yield differences between the sites were not differentiated by 
the accumulated stress degree days* As Implied earlier, the differences 
probably were due to Che effects of terrain slope and aspect* To examine 
this possibility data from the sites were segregated into four groups that 
represented good and poor drainage, and high, medium, and low levels of 
irradlance* These data are given in table 9. The slope and aspect for 
each site are given in table 1* Grouping the data as in table 9 shows 
that barley growth on south-facing slopes was generally better than on 
north-facing slopes, and that poorly drained soils can adversely affect 
plant growth* These data point out the complexities of the effect of 
slope and aspect on yield models that use remotely sensed data as inputs* 
Quantifying the effect of slope and aspect would not be feasible using 
only 1 year of data* Had Che experiment been conducted the previous year 
under conditions of drought, the south-facing slopes would probably have 
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been atreeeed end thus produce low yields » whereas the poorly drained 
sites may have been the only ones with adequate water, and thus would 
have had the highest yields* 

Comparison of Ground and Aircraft Acquired Surface Temperatures 

Ground-based temperature measurements offer the opportur ity to 
measure vegetation temperatures with a minimum of soil background 
Influence. This Is accomplished by viewing a target at an angle, after 
the plants have attained sufficient height. The crop water stress 
Index, and other temperature-based plant indices, depend on a measure of 
the temperature of only foliage. Scanners aboard aircraft and satelli- 
tes are essentially nadlr-looklng, and soil background may be 
observable. This problem was studied by Millard, Hatfield, and 
Goettleman (1979), and Millard et al., (1980). The latter compared 
ground-measured canopy temperatures with those obtained from an airborne 
scanner and concluded that, when 85Z or more of the soil surface Is 
covered by vegetation, airborne and ground measurements differed by less 
than 2*C. FurtheT- analysis Indicated that ground-based nadir measure- 
ments were nearly the same as the airborne data* 

When the soil surface Is moist. Its temperature will be closer to 
the vegetation temperature than when It Is dry* The frequent rain at 
the experimental site kept the soil surface moist much of the time. 

Toward the latter part of the season, plant senescence and higher sun 
elevations (more radiation striking the soil) caused the nadlr-looklng 
temperatures to exceed those measured at an angle (Figure 7). The effect 
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DAYS AFTER PLANTING 


Figure 7. The temperature difference between measurementa made holding the 
Infrared thermometer at an angle of about 30* from horizontal and 
when holding the instrument vertically, as a function of days 
after planting. 
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It nor* pronounced for tht afternoon aoeeuremente. After heeding (about 
day 120 to 130) an angular naaturamant vould be Influenced nora by haada 
chan by laavat, wharaaa the nadir view would Include aotcly laavat. Aa 
the leavaa aanaaca the nadir neaturaaenC %fould yield higher tamparaturaa 
than from an angle. Until 140 daya after planting, the angle-nadir tan- 
paraturaa differed by laaa than 1*C. Hatfield (1979) mtaeurad a maximum 
angle-nadir difference of -2*C whan planta covered 20 to SOX of the 
ground . 

An attempt waa made to ragiatar ground and aircraft data. The 
exact pixel in which the ground data were acquired could not be iden- 
tified in the airborne thermal imagery. The airborne data ware there- 
fore averaged over a 5 x 5 m pixel grid about the point where the ground 
aitea were thought to be located. Thia time conauming procedure waa 
done for seven flight daya. The reaulta are ahown in Figure 8. 
Regreoalon analyaia yielded a near zero intercept (-.34) and a alope of 
1.07, with an r^ of 0.97. 

A aecond compariaon waa made by averaging the aircraft pixel data 
over Che entire field and averaging the nadir-looking ground data for 
the 16 aitea. A total of 39 data aeta were available for thia purpoae. 
The reaulta are ahown in Figure 9. Again, a near zero (0.15) intercept 
waa obtained, aa well aa an identical alope (1.07) at the reglatered 
data in Figure 8. 

Theae compariaona lend confidence that methoda of detecting crop 
water atreaa and yield predictiona developed from ground data can be 
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Figure 8. Alrcrefu versus ground (nsdlr) measured temperatures. Alrcreft 
data were registered to the ground sites in time and space. 





CrvlG:ra>L F-Ti 
OF POOR QUALITY 


GROUND MEASURED TEMPERATURE 


igure 9. Aircraft versus ground (nsdlr) nsssured tsapsraturcs. Aircraft 

data wars averaged over the entire field. Ground data are averages 
of the 16 sites. 
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extended to aircraft-derived data. The cause of the slope being greater 
than 1 nay be due to Inadequate calibration of the hand-held radiometer 
or posbibly to an Inexact atmospheric correction being applied to the 
aircraft-derived temperatures. 

Spatial Integration of Aircraft-Measured Temperature 

In the previous section we established that a good correlation 
exists between ground- and aircraft-acquired temperature data. To 
relate aircraft to satellite data requires that the aircraft data be 
spatially Integrated. Data for the 1305 North-South flight of 16 Hay 
1978 was used for an example. The 2 x 2 m pixels were averaged into 
four ha square blocks (200 x 200 m) and shown in Figure 10 (top). The 
maximum, minimum, mean, and standard deviation are given for each block. 
Temperatures varied from a high of 44.3* to a low of 21.2“C. The mean 
temperature varied from 29.5* to 33.3*C. Further integration into 16 ha 
blocks (400 X 400 m) yielded a mean value range of 30.1* to 32.9*C. 
Integrating to 64 ha blocks (800 x 800 m) changed the range of the mean 
only slightly (30.3* to 32.7*C) (Figure 11). Integrating the data to 
256 ha (1600 x 1600 m), nearly the size of the field, yields a mean of 
31.5*C. The mean temperature of the largest block (256 ha) was within 
± 2*C from the highest and the lowest of the means of the 4 ha blocks. 

The spatial variation of temperature that occurs in a undulating 
field is rapidly attenuated as the resolution of the instrument 
Increases, as demonstrated in Figure 12. The S 3 rmbols in Figure 12 
represent the lowest and highest mean for each pixel size. The lines 
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Figure 10. Maximum, minimum, mean, and standard deviation for 
2m X 2m pixels averaged to 200m x 200m (4ha) pixels 
(top) and averaged to 400m x 400m (16ha) olxels 
(bottom). Data for aircraft flight at 130b, 16 May 78. 
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Figure 11. Maximum, minimum, mean, and standard deviation for 

2m X 2m pixels averaged to 80Qm x 800m (64 ha) pixels 
(top) and averaged to 1600m x 1600m (256ha) pixels 
(bottom). Data for aircraft flight at 1305, 16 May 78, 
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Figure 12. Mean temperature of pixels having the highest and lowest 

temperatures as a function of pixel size. Lines were drawn by 
Interpolating plots of the mean temperature versus the loga- 
rithm of pixel side length. 
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vere drawn by interpolating data from the plota of the mean pixel tem- 
perature and the logarithm of the pixel side length* Because of this 
interpolation, inferences about the data below 200 m side length are 
subject to question. With the length of the side of the resolution ele- 
ment of an HCMM pixel being 600 m, one can conclude that improving the 
resolution to 200 m would not significantly affect the estimation of the 
temperature within the field. Of course as the resolution improves 
below 200 m on a side, much more temperature detail becomes evident. 
Millard, Geottelman, and LeRoy (1981) reached the same conclusion from a 
different analysis of these data. 

Aircraft and Spacecraft Derived Imagery 

Figures 13, 14, and IS are images obtained from the low altitude 
aircraft data for the afternoon of 16 May 1978 and for the morning and 
afternoon of 20 May 1978. Imagery is available for other dates, but 
these were selected to correspond with HCMM data to be shown in sub- 
sequent figures. The three figures clearly show the temperature inhomo- 
genelty that occurs in undulating terrain. The northeast section of the 
field is the most rugged, and the Imagery shows that the temperature 
extremes occurred there. The morning image (Figure 14) shows a more 
uniform temperature distribution (note the expanded temperature scale). 
Temperature extremes differed by about 4*C. Figure 15 shows imagery for 
the afternoon of 20 May 1978. In general, surface temperatures were 
much warmer on 20 May than on 16 May. Temperatures were averaged over 
all pixels In the field for each Image and the 20 May temperature was 
7.S*C warmer than on 16 May. 
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Temperature Image of the Dunnigan agro-meteorologlcal 
experiment test field for 16 May 1978, in conjunction 
with a U-2 and HCMM overpass, approximate pixel size 
2m x 2m. 
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Figure 15. Temperature image of the field at 13J3, 20 May 1978, in conjunction 
with an HCMM overpass, approximate pixel size 2m x 2m. 
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In Figures 2, 13, and 15, the field boundaries In the upper right 
hand corner are quite distinct* The area outside these boundaries was 
not planted to a crop. Weeds and grasses were present because of the 
abundant rainfall. Comparing this area with the area Inside the field 
boundaries Indicates that the weeds and grasses were nearly the same 
temperature on 16 and 20 May. The large change of temperature within 
the field was due to ripening of the barley and the senescence of the 
leaves. This demonstrates that surface temperatures can be used to 
discriminate ripening. 

A U-2 derived Image for the afternoon of 16 May 78 Is shown as 
Figure 16* Temperatures within the field ranged from 30* to 34*C. This 
Is well trlthln the range of low altitude aircraft -derived temperatures 
(Figure 12). The U-2 pixel was about 70 m side length. The temperature 
range for this pixel supports the Interpolated data shown In Figure 12. 
Thus we can extend an earlier observation that, for fields of this typ*. 
(l.e., undulating terrain planted to the same crop), pixel sizes 70 m x 
70 m and larger yield essentially the same temperature range. 

HCMM Imagery at nearly the same time as the aircraft data are shown 
In Figures 17, 18, and 19. Location of the site was extremely difficult 
and we are not completely confident that the nine pixels used were 
totally within the field* Another complication was the calibration of 
the HCMM radiometer. Barnes and Price (1980) reported the pre-launch 
calibration, and also a post-launch calibration over a lake In New 
Mexico. The post-launch values were about 5.5*C below the pre-launch 
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perature image of the field at 1337, 16 May 1978, 
aincd with the U-2, approximate pixel size 70m x 70m. 
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Figure 18. Temperature image obtained with the HCMM satellite at 

approximately 0230, 20 May 1978, approximate pixel size 
600m X 600m. 
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v«lu«s. Figures 17-19 srs In terms of the poet-leunch celibretlo;'.. We 
obtslned tenpereture velues from where we believe the field '«as located. 
These date are compered with the aircraft data In Figure 20. The poet- 
launch data fall well below the elrcraft late. When 5.5*C is edded to 
the HCMM values, they fall much closer to th&^e from the aircraft. 

Three data points are not sufficient to calibrate the HCMM date with the 
elrcraft data hut they are sufficient to show chat the pre-launch 
calibration was superior to the post-leunch one. It la possible that 
the atmospheric correction epplled to the HCMM data was coo high at 
night and too low during the day. The U-2 veluc falls nearly on the 1:1 
line. 

CONCLUSIONS 

A crop water stress Index was developed from energy balance con- 
siderations to quantify the effect of stress on plants by using remotely 
sensed canopy temperatures. 

Aircraft derived surface temperatures are well correlated with 
ground based temperature ia.'^eurement8. Thus, ground-based measurements 
can be economically used to develop strategies for soil moisture and 
crop water stress evaluations that can be readily adapted to aircraft 
derived data. 

Integrating aircraft derived data from 2 a x 2 a pixels to 200 a u 
200 a pixels reaoved spatial variation of temperature to the point that 
further Integration to 800 m x 800 m pixels showed essentially the same 
temperature detail. U-2 data indicate that pixels greeter than 70 m x 
70 m yield results similar to the 200 m x 200 m pixel. These results 
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AIRCRAFT MEASURED TEMPERATURE 


Figure 20. Satellite versus aircraft measured surface temperatures. The 
square symbol Indicates the U-2 derived temperatures. The x’s 
and circles represent HCMM data using the pre- and post>launch 
calibrations, respectively. 
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suggest that there is little advantage to decreasing pixel size from 800 
m X 800 m to 70 o x 70 m as far as assessing temperatures over undulating 
terrain is concerned. The smaller pixels, however, would make it much 
easier to locate specific areas only 9 large pixels in size. 

The pre-launch calibration of the HCMM radiometer yielded tempera- 
tures closer to aircraft derived temperatures than did the post-launch 
calibration. Using the pre-launch calibration, temperatures were higher 
at night and lower during the day than were the aircraft temperatures. 
This raises a question as to the calibration of the HCMM radiometer or 
to the adequacy of regionally acquired radiosonde data used in conjunc- 
tion with algorithms to adjust for atmospheric effects. 


Intensive ground and aircraft experiments that are geared to the 
growth and development of plants should not be committed in conjunction 
with a satellite program until the satellite is launched and is 


operational. All biological field experiments are subject to the 
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